arXiv: 1509.04741 v2 [cond-mat.mes-hall] 7 Dec 2015 


RKKY interaction of magnetic impurities in Dirac and Weyl semimetals 

Hao-Ran Chang/ Jianhui Zhou/’* Shi-Xiong Wang/ Wen-Yu Shan/ and Di Xiao^’^ 

^Department of Physics and Institute of Solid State Physics, 

Sichuan Normal University, Chengdu, Sichuan 610066, China 
^Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA 

(Dated: December 8, 2015) 

We theoretically study the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between mag¬ 
netic impurities in both Dirac and Weyl semimetals (SMs). We find that the internode pro¬ 
cess, as well as the unique three-dimensional spin-momentum locking, has significant influences 
on the RKKY interaction, resulting in both a Heisenberg and an Ising term, and an additional 
Dzyaloshinsky-Moriya term if the inversion symmetry is absent. These interactions can lead to 
rich spin textures and possible ferromagnetism in Dirac and time-reversal symmetry-invariant Weyl 
SMs. The effect of anisotropic Dirac and Weyl nodes on the RKKY interaction is also discussed. 
Our results provide an alternative scheme to engineer topological SMs and shed new light on the 
application of Dirac and Weyl SMs in spintronics. 
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Three-dimensional (3D) Dirac semimetals (SMs) [1] 
are topological states of matter and can be seen as bulk 
analog of graphene. Their conduction and valence bands 
with linear dispersion touch each other at a finite num¬ 
ber of points, called the Dirac nodes, in the 3D Brillouin 
zone. Dirac nodes are fourfold degenerate, protected by 
both time-reversal symmetry (TRS) and inversion sym¬ 
metry. Breaking either symmetry in Dirac SMs leads 
to Weyl SMs [2], which host Weyl nodes. These Weyl 
nodes can be viewed as effective magnetic monopoles in 
the momentum space [3], acting as the source and drain 
of the Berry curvature field [4]. This nontrivial topol¬ 
ogy can lead to exotic superfluid [5] and superconduct¬ 
ing phases [6], unique Fermi arc state [7], helical spin or¬ 
der [8], various novel electromagnetic responses, such as 
the chiral anomaly [9, 10], the chiral magnetic effect [11], 
negative magnetoresistance [12], and the chiral Hall ef¬ 
fect [13]. 

Many aforementioned exotic phenomena rely on the 
separation of Weyl nodes in momentum space due to 
the intrinsic TRS breaking. So far, angle-resolved pho¬ 
toemission spectroscopy and magnetotransport measure¬ 
ments have identified (Bii_:jIn 2 ;) 2 Se 3 [14], NaaBi [15], 
and Cd 3 As 2 [16] as Dirac SMs, and noncentrosymmet- 
ric transition-metal monosphides TaAs, NbAs, NbP and 
TaP as Weyl SMs [17-22]. There are, however, few ex¬ 
perimental realizations of TRS breaking Weyl SMs with¬ 
out Landau quantization [23]. Even though a magnetic 
field can break the TRS, it inevitably couples to both 
the spin and orbital motion of electrons through Zeeman 
splitting and Landau quantization, respectively. On the 
other hand, magnetic doping technique has recently been 
utilized in experimental implementation of the quantum 
anomalous Hall effect in thin films of topological insu¬ 
lators [24]. Thus, one may naturally wonder whether 
Weyl SMs without TRS can emerge from Dirac and time- 
reversal invariant Weyl SMs through introducing mag¬ 


netic dopants rather than applying a magnetic field. 

Motivated by the above observation, in this Rapid 
Communication we study the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction between magnetic dopants 
in both Dirac and Weyl SMs. We find that the internode 
process, as well as the unique 3D spin-momentum lock¬ 
ing, has significant influences on the RKKY interaction, 
resulting in both a Heisenberg and an Ising term, and 
an additional Dzyaloshinsky-Moriya term if the inver¬ 
sion symmetry is absent. These interactions can lead to 
rich spin textures and possibly ferromagnetism in Dirac 
and TRS-invariant Weyl SMs. The effect of anisotropic 
Dirac and Weyl nodes on the RKKY interaction is also 
discussed. Our results provide an alternative scheme to 
engineer topological SMs and shed new light on the ap¬ 
plication of Dirac and Weyl SMs in spintronics. 

In general, a pair of Weyl nodes of opposite chirality 
can be described by the following Hamiltonian 

Hg = x[vF(T ■ (k - xQ) P ooQo] , (1) 

where k is the wave vector, vp is the Fermi velocity, 
and X = refers to the chirality of the Weyl nodes. 
If {Q,Qq) = 0, the two Weyl nodes overlap with each 
other and the Hamiltonian Hg describes Dirac SMs [see 
Fig. 1(a)]. For noncentrosymmetric Weyl SMs that pre¬ 
serve the TRS, Q must be zero and Qg can be nonzero. 
As a result, the two Weyl nodes are located at the 
same k point but can have different energies as shown 
in Fig. 1(b). Hence for a given carrier density, there are 
two unequal Fermi wave vectors. On the other hand, for 
Weyl SMs with broken TRS but with inversion symme¬ 
try, we have Qg = 0 and Q yf 0; the two Weyl modes have 
the same energy but reside at different k points ±Q in 
the Brillouin zone [see Fig. 1(c)]. Here the Pauli matrices 
O' = [<Jx,cry,az) refer to the real spin degree of freedom 
of electrons. They may also refer to pseudospin degree 
of freedom, which will be discussed later. 
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FIG. 1. (Color online) Low energy spectra for Dirac SMs (a), 
noncentrosymmetric Weyl SMs (b), and Weyl SMs without 
TRS (c), respectively. The color of cone indicates the chirality 
of Weyl nodes, while the blue plane is the Fermi level £f = 0. 


We assume that the interaction between the 3D itin¬ 
erant Weyl ferirrions and a irragnetic iirrpurity Si located 
at R^ can be expressed as the standard s-d interaction 
Hamiltoniair Hj = (Jtq + XTx)Si ■ crSir — Ri), where J 
and A refer to the strength of the s-d exchange interac¬ 
tion in the intranode process and the internode process, 
respectively. The identity matrix tq and Pauli matrix 
act on the chirality space. 

Accordiirg to the second order perturbation the¬ 
ory [25], at zero temperature the RKKY interaction 
between two magnetic impurities mediated by Weyl 
fermions is given by H^kky = de Tr[( Jtq -b 

Xtx)Si ■ (TQ{R-,e){jTQ + Xtx)S 2 ■ a-g^-R-e)] with R = 
R 2 — i?i, where £f is the Fermi energy and Tr means 
a trace over the spin and pseudospin degree of freedom 
of itinerant Weyl fermion. The Green’s function in the 
energy-coordinate representation is given as g{R;e) = 
G+{R-,e) 0 G-{R\e). After some algebra, we find the 
RKKY interaction [26] 

d^RKKY= E [J^S^^,+X^{1-S^^,)]S?S^ 

a,P,X,X' 

xlm|- J deTT:[aaG^{R]£)apG^'{—R]£)\^ , 

( 2 ) 


which includes contributions from both the intranode 
process and the internode process. 

We now apply the above formalism to study the RKKY 
interaction in Dirac and Weyl SMs. We note that close 
to the Weyl nodes, the effective Hamiltonian is propor¬ 
tional to fe • (T, which can be viewed as the 3D counter¬ 
part of graphene in real spin space. It is in contrast with 
the surface state of 3D topological insulators (TIs) [27- 
29] , in which spin and velocity are perpendicular to each 
other. As we show below, this kind of hedgehog spin tex¬ 
ture around Weyl nodes (spin is aligned with momentum) 
can have significant effects on the magnetism of magnetic 
impurities. 

Let us first consider the isotropic Dirac SMs, described 
by Huik) = x^pk ■ cr with the energy dispersion = 


-iiVFk. The effect of anisotropic energy dispersion will be 
discussed later. The Green’s function corresponding to 
H'£i[k) in momentum space takes the form G“^(fe;e) = 
(e + ir])ao — H]j{k), where ry is a positive infinitesimal. 
G^{±R-,£) can be obtained from G^{k;£) through in¬ 
tegrating over the momentum near the Weyl node y, 
G^{±R;£) = / (f^Gx(fc;e)exp(±ife • J?). Garrying 
out the integration over k leads to (more details are pre¬ 
sented in the Supplemental Material [30]) G-^{±R;£) = 
aoGo{R-,£) ± xcrjGi^(i?;e), where the Green’s functions 
are defined as Go(i?; £) = ^ exp)!^^] and Gh{R; e) = 

iTrvFii^ (1 — z^e) exp[f^£], with the dimensionless parame¬ 
ter = £R/vf- In the above calculation, we take R to 
be aligned to the j axis, i.e., R = Rsj. 

After lengthy but straightforward calculations, we ob¬ 
tain the RKKY interaction for the Dirac SMs (the de¬ 
tailed calculation can be found in the Supplemental Ma¬ 
terial [30]) 

Hrkky = fS{R, ^f)Si ■ S 2 + FD„g(R, ^F)SiSi , (3) 

where the range functions for the Heisenberg and Ising 
terms are given by 

= - {-/^[(S - 2^p) cos(2^f) + sin(2^F)] 

—2A^[cos(2^f) + sin(2^F)]} /Stt^vfR^ 

(4) 

Cf) = -{J^ - x^) [(2^^ - 5) cos(2eF) 

—6^_f sin(2^F)], (5) 


with ^F = £fR/vf- Here the superscript D stands for 
Dirac SMs. One can clearly see that if = A^, the 
contribution to the Ising term from the internode pro¬ 
cess cancels that from the intranode process exactly, i.e., 
Kj^ing = 0. This cancellation is due to the restoration of 
the spin rotation symmetry by the internode process. 

At long range {^f 1) and for finite Ef, the RKKY 

interaction reduces to a simple form 




D 

RKKY 


J^£|,cos(2^f) 


( 6 ) 


which reflects the long-range and oscillatory nature of the 
RKKY interaction. The ferromagnetism (FM) or antifer¬ 
romagnetism (AFM) of magnetic impurities depends on 
both the concentration of impurities and the carrier den¬ 
sity of Dirac fermions through R and £f, respectively. 

Interestingly, for the intrinsic case (ef = 0), the 
RKKY interaction becomes nonoscillatory: 


rrD 

^RKKY 


-(3 - 2A^) 

Stt^vfR^ 


Si- 82 . 


(7) 


We can see that if = A^, the exchange coupling is al¬ 
ways ferromagnetic. For sufficiently large impurity den¬ 
sity, this will lead to a spontaneous magnetization, which 
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FIG. 2. (Color online) The exact range functions of RKKY in¬ 
teraction of the noncentrosymmetric Weyl SMs (i^^^/C) as a 
function of the reduced parameter C) at the Fermi energy ef = 
15Qo (a) and ef = O.lQo (b), with C = 

The inset in each panel for shows an evident beat¬ 

ing feature for each term in this RKKY interaction. 

then drives Dirac SMs into Weyl SMs with broken TRS. 
Consequently, a nonzero anomalous Hall conductivity 
proportional to the separation of the two Weyl nodes 
in momentum space will appear [31]. The ferromag¬ 
netic transition temperature Tc depends on the specifics 
of both magnetic impurities and host materials, which 
need detailed first-principles studies [32]. It is clear that 
its spatial dependence as 1/R^ differs from 1 /R^ for the 
intrinsic graphene [33]. 

Next we consider noncentrosymmetric Weyl SMs, in 
which the two Weyl nodes with opposite chirality have 
different energies {Qq ^ 0). For a given carrier density, 
the magnitudes of Fermi wave vectors of these two Weyl 
nodes are distinct, therefore ^ with = {ep — 
xQo)R/vf- Following the same procedure, we obtain the 
RKKY interaction for the noncentrosymmetric Weyl SMs 


pressions of these range functions Ti™g, ^"dm) 

given in the Supplemental Material [30]. 

Figure 2 shows the range functions to two sets of pa¬ 
rameters. As we can see, these range functions display 
a damped oscillatory behavior with increasing distance 
R, with each term dominating in different regimes of the 
parameters e p and R. In addition, these range functions 
oscillate with two distinct periods and form a beating 
pattern. The beating feature, originated from the two 
unequal Fermi wave vectors, manifests itself by multiply¬ 
ing each of the range functions in Eq. (8) by as shown 
in the insets of Figs. 2(a) and 2(b). The measurement 
of beating period can be used to determine the energy 
difference Qq [35]. This beating structure does not occur 
in the typical surface state of TIs where there is only one 
Fermi circle [27-29]. 

The DM term is essential for realizing spiral spin states 
and skyrmions, and also has potential applications in 
spintronics. It may also provide some hint to understand 
the recent experimental observation of the TRS-breaking 
Weyl metal in YbMnBi 2 [23]. In this material, although 
there is a global inversion symmetry, in each layer the 
inversion symmetry is broken, which could give rise to 
a nonvanishing DM term. The combination of the DM 
interaction and the AFM order of Mn can cause a canted 
AFM order with a nonvanishing net magnetization ob¬ 
served in the experiment [23]. 

Finally, we discuss TRS-breaking Weyl SMs, in which 
the Weyl nodes with opposite chirality reside at different 
k points in momentum space ±Q. We further assume the 
inversion symmetry remains intact, thus the DM term 
does not appear. We take into account the effect of the 
separation of Weyl nodes in the internode process. The 
corresponding RKKY interaction can be obtained from 
that of Dirac SMs in Eqs. (4) and (5) by replacing with 
A^ cos(2Q ■ R). It can be seen that the internode process 
gives rise to an oscillating term proportional to ^ cos(2(5- 
R), which is absent in the typical surface state of TIs [27, 
28]. Because of the long-range and oscillatory nature of 
the RKKY interaction, for a large momentum separation 
2Q, the part from the internode process will vanish after 
averaging over position and thus does not contribute to 
the net magnetization. Therefore, the intranode process 
dominates the RKKY interaction in this case. At long 
distance, we have the approximate RKKY interaction 


i^RKKY = OSl ■ S2 + C) 

X SiSi + C)(‘^i X S2)j , (8) 




RKKY 


j'^Ep COs(2^i7) 
4:Tr^VpR^ 


(Si ■ S2 - SiSi) , 


(9) 


with C = QqR/vp and the jth component of the spin 
is along the direction connecting the two impurities. 
Here the superscript TR stands for the TRS invariant 
Weyl SMs. We can see that the RKKY interaction con¬ 
sists of three terms, namely, a Heisenberg term, an Ising 
term, and a Dzyaloshinsky-Moriya (DM) term [34]. This 
is similar to the surface state of TIs [28] . The specific ex- 


where the superscript I refers to the centrosymmetric 
Weyl SMs. Since the asymptotic range functions for the 
Heisenberg and Ising terms exhibit the same oscillatory 
behavior but with a phase differing by tt, the Ising term 
always cancels the jth component of the Heisenberg term, 
leading to the AY-like spin model for j = z. On the 
other hand, for the intrinsic case the RKKY interaction 
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becomes 

which could realize various spin models, such as the 
XXZ-like spin model for j = z. The resulting spin con¬ 
figurations of impurities {XY- and XXZ-\\kB spin mod¬ 
els) can be accessed by a variety of experimental tech¬ 
niques, such as neutron scattering technique. 

In reality, almost all experimental realizations of Weyl 
or Dirac SMs (and theoretically conjectured systems, 
such as [31, 36-38]) possess a strongly anisotropic single¬ 
particle dispersion. Here we consider the effect of 
anisotropy on the RKKY interaction. The effective 
Hamiltonian for fermions near the anisotropic Weyl node 
X takes the form Hq = xi'^xkx^x + Vykyay -I- v^k^tTz), 
where the Fermi velocities Va in different direction are 
different, \vx\ ^ 1%] ^ \vz\- To simplify our discussion, 
we set Vj > 0 with j = x,y,z. The corresponding en¬ 
ergy spectrum is given as = ±{vlk'^ + Vyky + vlkiy^'^. 
It is instructive to make the following transformation as 
(viQi^Viki^vpRi) = {vFQi,VFki,ViRi), such that we re¬ 
late the Green’s function of an anisotropic Weyl node 
in the energy-coordinate representation to that of an 
isotropic one as G'^(±i?;e) = XaG^{±R-,£), with = 
Vp/ivxVyVz)- We note that the above transformation 
converts a Fermi elliptic sphere into a Fermi sphere but 
preserves the volume of Fermi sphere. This allows us to 
connect the RKKY interaction of anisotropic Dirac SMs 
to the counterpart of the isotropic ones 

HjiKKY{R) = X\HjiKKY{vFR/Vi) , (11) 

which implies that the anisotropy of Dirac or Weyl nodes 
must lead to an anisotropic RKKY interaction. It is clear 
that for Ua; = Vy = Vz (A^i = 1 ), the above expression in 
Eq. (11) can reduce to that of the isotropic one. 

Before drawing conclusions, we briefly discuss the 
pseudospin case, in which the Pauli matrices in Eq. (1) 
may refer to pseudospin degree of freedom such as orbital 
index. This pseudospin case is similar to graphene. In 
the absence of the spin-momentum locking, the RKKY 
interaction only contains the conventional Heisenberg 
term [39]. Hence, the RKKY interaction directly asso¬ 
ciates with the Fourier transform of the static density- 
density response function [40]. We leave all the specific 
expressions of RKKY interaction to the Supplemental 
Material [30]. It should be emphasized that the RKKY 
interaction in this intrinsic Dirac SMs is also a nonoscil- 
latory Heisenberg term Hb.kky = 4)^3 Si ' S 2 , which 
allows the spontaneous magnetization of magnetic impu¬ 
rities for jA] > (2/3)^/^] J]. Compared with Eq. (7), the 
factor of 2 comes from the degeneracy of real spin degree 
of freedom. 

In summary, we have studied the RKKY interaction 
between magnetic impurities in Dirac and Weyl SMs. We 


found that it is possible to realize a spontaneous magne¬ 
tization in these systems. The RKKY interaction in gen¬ 
eral contains the Heisenberg, Ising, and DM terms, which 
can give rise to rich spin textures of the impurities. These 
findings provide an alternative scheme to engineer topo¬ 
logical SMs and pave the way for the application of Dirac 
and Weyl SMs for spintronics. 

Note added. Recently, a related paper appeared, in 
which some results on the anisotropic Weyl SMs have 
also been obtained [41]. Also recently, the Kondo effect 
of a single magnetic impurity in the Weyl SMs has also 
been discussed in [42, 43]. 
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SUPPLEMENTARY MATERIAL FOR”RKKY INTERACTION OF MAGNETIC IMPURITIES IN DIRAC 

AND WEYL SEMIMETALS” 

This supplemental material is organized as follows. In Sec. I, we present the general expression of RKKY interaction 
in terms of the Green function. In Sec. II, the Green function in the energy-coordinate representation is calculated 
detailedly. In Sec. Ill, we further calculate the RKKY interactions for 3D Dirac SMs, Weyl SMs without time reversal 
symmetry (TRS), and noncentrosymmetric Weyl SMs in both cases with and without spin-momentum locking. In 
Sec. IV, we investigate the impact of anisotropy of Weyl or Dirac nodes on RKKY interaction. In Sec. V, we give 
some useful integrals for calculating the range functions. 


FORMALISM OF RKKY INTERACTION 


Following the second order perturbation theory [25] , we express the general RKKY interaction between two magnetic 
impurities mediated by itinerant fermions in terms of the Green function as 


Hkkky = -Im 

TT 


de Tr 


(Jto -I- Xtx) Si ■ a-g (R, er) (Jtq -I- Xt^) 82 - crQ (-R, e) 


( 12 ) 


where ep is the Fermi energy, Tr means a trace over the (pseudo)spin degree of freedom of itinerant Weyl fermion, the 
identity matrix tq and Pauli matrix Tx act upon the chirality space, J and A refer to the strength of the s-d exchange 
interaction in the intranode process and the internode process, respectively. The Green function in energy-coordinate 
representation is given by 


g (±r,£) 


/G+(±R,e) 0 \ 

0 G-{±R,e)J^ 


(13) 


where (±R,£) is the Fourier transform of the Green function in momentum space G^(k, e). Inserting Eq. (13) into 
Eq. (12), we obtain the RKKY interaction 


Rrkky = E + A^(l - '5xx')] {R; sp), 

x,x' ,a,l3 

where the functions in different channel are given by 


T^^iR;ep) = - 


1 


d£Tr [o-aGx (R,£) (-R,£)] ■ 


(14) 


(15) 


It should be noted that Eq. (14) is our main expression of RKKY interaction, and the evaluation of {R;ep) is 
left to the next part. 


DERIVATION OF THE GREEN FUNCTIONS 

We at first calculate the Green function of the Dirac SMs in energy-coordinate representation, which is the Fourier 
transform of that in momentum space G~^ (k, e) = (e -|- ir]) uo — X'^pcr ■ k. 


G^ (±R,£) — 




G^ (k, e) exp(±ik ■ R) = 


d^k eao iyuFcr • k 

'o - —pi - 5 - exp (zk • R) , 

(27r)3 (e + ir])^-el 


(16) 


where Sk = ivpk. In general, the momentum k can be decomposed into the directions parallel and perpendicular to 
e/{ as k = k|| -|- kj_ = (k • ep)eii + {ep x k) x ep with ep = R/|R|. Because of kp- ep = 0, the integration over 
angular variables causes the part of kj_ in G^^- (R; £) to vanish, we then have 

Gx(±R, £) = croGo(i?; £) ± ycr • epGp{R-, e). (17) 

Therefore we only need to evaluate two integrals, Go{R;e) and Gp{R;e). The former is calculated as follows, 

d^k exTp{ikRcos6) e k^dk 


Go(R; e) — e J 


{27r)^ {e + irj)^ — Vpk'^ Jq k'^ — {e + irj)'^/Vp 


sinO exp{ikR cos 0)d6, (18) 
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where 9 is the angle between vector k and vector R. Integrating over 6 leads to 

k [exp (ikR) — exp (—zfci?)] dk 

VpP I 




— (e + ir])'^/v1 


Due to the symmetry of the integrand of k, the integration of Gq{R;£) can be recast as 


Go{R; e) — —- 


fcexp (ikR) dk 


iAir^VpR — (e + iri)'^/vp 

After carrying out the contour integration, we arrive at the result, 

Go{R;£) = - exp (ieR/vp) = = G'ofe), 

where = eR/vp is a dimensionless quantity. Since the latter one Gp(R;£) satisfies the relation 


Gr{R; e) = vp 


k • br exp(ik • R) 


(27r)3 (e + ir])^ — Vpk^ 
we immediately get the latter 


f d^k k cos 9 exp (ikR cos 9) ■'^^d 

J (2tt)^ (e + — Vpk"^ * e dR ° ’ 


Gr(R]£) = -- — ^—^i{l-i£R/vp)exp{i£R/vp) = --—(1 - exp =GR(^e). 


AttvpR'^ 


iTrypR"^ 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Note that the procedure above is also applicable for the Weyl SMs. Taking into account the separation of Weyl nodes 
in the internode process, we get the corresponding Green function for Weyl SMs without TRS in energy-coordinate 
representation 


/ d^k f d^k 

—^Gx(k,e)exp(±i(k-l-xQ)-R) =exp(±ixQ-R) / —(k, e) exp(±ik • R), (23) 

^ZTrj J (^ZTrj 

which is the counterpart of the Dirac SMs multiplied by a factor of exp (±ixQ • R)- For noncentrosymmetric Weyl 
SMs, the corresponding Green function in the energy-coordinate representation can be directly obtained from that of 
Dirac SMs by replacing e with e — xQo- 


RANGE FUNCTIONS AND THEIR ASYMPTOTIC BEHAVIORS 


To simplify the presentation, we only show the detailed calculation of the RKKY interaction of isotropic Dirac SMs. 
Substituting Eq. (17) into Eq. (15), we get the corresponding RKKY interaction 


HRKKviR-.eF) = Im ^ V + 

X.X'.aA 


with the functions 


Rap (R'^^f) = J aa(aoGo(R;£) + X(^3GR{R;£)^(7fj(aoGo(R]£) - x'(^jGR(R;£)^ 


(24) 


(25) 


where we have assumed br = Bj , which implies cr • br = Uj . For the remainder of this section, we compute the RKKY 
interactions for the cases with and without spin-momentum locking, respectively. 


Case with spin-momentum locking 


For the case with spin-momentum locking, the trace term in Eq. (25) can be calculated as follows. 


Tr 


CTa(o-oGo(R;e) + xo'jG'ij(i?;e)^cr;g(^(ToGo(R;e) - x'<^jGRiR;£)^ 


=Go(i?; e)Go(i?; £:)Tr (CTQcroa-/3cro) - Xx'Gr(R-, £)Gr(R-, £:)Tr (aaCCjapaj) 

-I- xGr(R-, £)Go(R] e)Tr (aaCCjcrpao) - x'Go(R; £)Gr(R; e)Tr (aaaoapaj). 


( 26 ) 



















Substituting it into Eq. (25), we have the functions 


ttxx 

•^a/3 


^R-ep) = -i (^) 1^" [2S^p(Go{OGo{0 + Xx'Gr{OGr{0) 

- iS^jSp.xx'GRiOGRiO - 2*£a/3,(x + x')Go{OGr{0 


with = EfR/vf- Making use of Eq. (24) and summing over x, X^ ^-^d /?, we get the RKKY interaction of the 
Dirac SMs 


^fRKKY(^;^^’) = ■ S2 + FE,^^{(F)SiSi, 

where the range functions are given as 

r^F 


(27) 


IT’D 


((f) = -i (^) ImJ' d(l [j2(g„(0G„(0 + G*(0Gk(0) + A"(g„(0G„(0 - Gk(0Gk(0) 


f\r) 


1 


AfvfR^ 

rip 


4 Im 


J2 (2l2i^F, 2) + 2iF{CF,2) - Io(Cf, 2)) + (^-2iFiCF,2) +Ioi^F, 2)) 


FEinsi^F) = -I (^) Im f\^8{X^- J2) Gfi(CGR(C 

TV \ IX / J _oo 

- _i 

TT V i? / 


j 8 (A^ - J 2 ) Im[X 2 (^F, 2 ) + 2 ili(^f, 2 ) - Xo(ef, 2 ) 
where we have utilized Eqs. (20) and (22). Making use of the basic integral I„(^F,a) given in Sec. , we finally obtain 
FE{^f) = — ^(3 — 2 ^p) cos( 2 ^iA) + 4^iA sin( 2 ^iA)^ — 2 A^ ^cos( 2 ^iA) + sin( 2 ^ir)^ Istt^vfR^, (28) 

FEsingi^F) = A') [(2^1. - 5) cos(2eF) - e^F sin(2eF)] /Stt^fR^. (29) 

Note that the Ising term FEi^^i^F) vanishes for homogeneous coupling case A = J. The asymptotic behavior of range 

function (^f) at long range (^f 2 > 1) reduces to (^f) ~ J^Ef oos(2^f) j‘^ f^vfR^ = cos(2fF) j 

The procedure above is also applicable for the Weyl SMs. For the Weyl SMs without TRS, (E; Ef) in Eq. (25) 

is substituted by exp[i(x — xOQ ' {R',Ef)- After summing over x, X^ ct, and /3, the corresponding RKKY 

interaction becomes 

AIrkky iR-,EF) = F^i^F)S, • S 2 + E/,i„J^F)SlS^2> 
where the range functions are given by 

F^i^F) = — <^^^(3 — 2^p) cos(2^f) + d^F sin(2^F)^ — 2A^ cos (2Q ■ R) ^cos(2^f) + Cf sin(2^F)^ 

Fiingi^F) = A" cos (2Q • R)) [( 2^2 _ 5 ) cos(2eF) - G^f sin(2^F)] / 8tt\fRE 


(30) 


Stt^vfR^ 


(31) 

(32) 


Two remarks on the range functions are in order here. First, F^{^f) and F^ngi^F) can be obtained from Eg (^f) 
and Ej°j^g(^F) by replacing A^ by A^ cos (2Q • R). Second, the internode process gives rise to an additional oscillating 
term proportional to cos (2Q • R). Because of its long range and oscillatory nature of the RKKY interaction, for a 
large momentum separation 2Q, the part from the internode process will vanish after averaging over position. Hence, 
the intranode process dominates the RKKY interaction in the Weyl SMs without TRS, leading to a non-vanishing 
Ising term. Therefore the range functions can be expressed as 


F^i^F) = - 2^f) cos(2,c^) -I- 4^f sin(2,c^)^ j8 FvfR^, 

Fisingi^F) = {2 Ef - 5) cos(2,c^) - 6 ^f sin(2^F) j8 tF‘vfRE■ 


(33) 

(34) 
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The asymptotic range functions at long range ^ 1) become 

-Ph(Cf) ~ COs{2^f) j'^TT^VpR^ = J^£^f cos(2,fi;’) j 

^Ising(CF) « -J^^F COs{2^f) jA. f'^VfR^ = -J^Ep C0s(2^f) j^' k'^V%R^. 

For the noncentrosymmetric Weyl SMs, {R-ep) in Eq. (25) turns out to be 




1 

TT 



de Tr 


^Oi 


X(70 


(aoGo{R] £ - xQo) + xf^jGuiR] £ - xQo)) 
{aoGo{R-, £ - x'Qo) - x'<^jGR{R; e - x'Qo)) 


After summing over the indices for spin and chirality, the corresponding RKKY interaction becomes 
^R^KY(i?; ep) = FS^i^FX)S^ • 52 + OSjSi + fS^Xf, () (5i x 52),-, 

where = QqR/vf and the three range functions are given as 


(35) 

(36) 


(37) 

(38) 


(39) 


Fu'^i^FX) 

^’iIi^g(eF,C) 

FS^i^FX) 


^(3 — 2^p — 2(^) cos(2^i^’) cos (2() + 4^f( sin(2^i7’) sin (2^) + 4^p sin(2^i7’) cos (2() 

+ 4^cos(2^i^’) sin (2(^)^ — 2A^^cos(2^i7’) + C_F sin(2^i7’)^ j^F^ypR^, (40) 

(^(2^f + 2(X' — 5) cos(2^i;’) cos (2^) — 4,^^^ sin(2^i7’) sin (2Q — Q^p sin(2^F) cos (2^) 

— 6C cos(2^i;’) sin (2(^)^ + ^(5 — 2^^ + 2(X) cos{2^f) + 6^f sin(2^i7’)^ jSn^ypR^, (41) 

(^(2^f + 2(X — 2) cos(2^i^’) sin (2^) + 4^f( sin(2^i;’) cos (2^) — 4^^ sin(2^i7’) sin (2^) 

+ 4()cos(2,c^) cos (2())^ — 2 A^Ccos(2^f) js^^ypR^. (42) 


It is clear that for C = 0, i.e., the Dirac SMs case, the DM term ^'^^(^iT’.O) vanishes, and the range functions F™(^i;’^C) 
and F^^^XfX) reduce to F^Xp) and F^^^^Xp), respectively. In addition, the range functions of RKKY interaction 
for the noncentrosymmetric Weyl SMs oscillate with two distinct periods and form a beating pattern. Note that 
these behaviors can be roughly described by the asymptotic expressions at long range. For ^p :s> ( and ^p ^ 1, the 
asymptotic range functions become 

FS^XfX) - JHI cos(2ej’) cos (20 I4fHfRX 
F^I^XfX) « -Cf cos(2eF) (J^ cos (20 - A^) /4p\fRX 
^dm(CdC) ~ -J'^ip cos(2^f) sin (20 j4 p'^vfR^. 


For ( ^ ^p and C ^ 1; the range functions read 

F^^XfX) « cos(2Cf) cos (20 /4f\fRX 

F^sE^^XfX) « -C cos(2eF) (J' cos (20 + A^) /47t\pRX 
FSmi^FX) ~ -J^C^ cos(2,Cf) sin (20 j4 p'^vfR^. 


It should be emphasized that these two different asymptotic range functions account for the distinct beating pattens 
in Fig. 2(a) and 2(b) in the main text, respectively. 










10 


Case without spin-momentum locking 

For the case without spin-momentum locking, the trace term in {R^ep) reads 
, (toGo{R; e) + xtjGr{R; e)^a 0 (toGo{R; e) - x'fjGR{R; e)) 


Tr 


=GoiR; e)Go{R-, e)Tr(craCr;3)Tr(foTo) - Xx'Gr{R-, e)Gfi(i?; e)Tr(cr„cr/3)Tr(fjfj) 

-I- xGr{R-, e)Go{R; e)Ti{aaCrp)Tr{TjTo) - x'Go{R-,e)GR{R\e)TT{aaap)Tr{ToTj), 


(43) 


Where the identity matrix tq and Pauli matrix fi here refer to other pseudospin degree of freedom rather than the 
chirality. Substituting the above trace term into the definition of {R; er), we have 

UR-.sf) = -I (^) d^iS^p{Go{OGo{0-Xx'GR{OGR{0)- (44) 




After summing over the indices for the spin and chirality degree of freedoms, we get the corresponding RKKY 
interaction 


'^rkky(4?; £^f) — FEi^F)Si ■ S 2 , 


(45) 


where the tilde is used to distinguish the range function of the case without spin-momentum locking from that with 
spin-momentum locking. Utilizing the explicit expressions of basic integral In{^F,ci) evaluated in Sec. , we finally 
obtain the range function, 


Fh'(Cf) = 4J^ ^cos(2^f)-|- sin(2^F)^ — ^2(3 — 2^1^) cos(2^i^’)-I-sin(2^i^’)^ j^ f^vrR^ 

Similarly, the RKKY interaction for the TRS-breaking Weyl SMs reads 

'^rkky(^Uf) = -Fh(Cf)‘S'i ■ S2, 

where the range function is 

^h(Cf) = 4J^^cos(2^i?)-I-sin(2^F)) -A^cos(2Q-R) ^2(3 - 2^|) cos(2^i?)-|-8 ^f sin(2^i?)^ 
and that for the noncentrosymmetric Weyl SMs, 

Hll^^{R-eF) = Fl^{^F)S,-S 2 . 

where the range function has the form 


(46) 

(47) 

^8F3^;Fi?^ (48) 
(49) 




-TR 


(^f) = 4 ^cos(2,^f) cos(2C) -I- Cf sm(2,^F) cos(2() -|- C , cos(2,^f) sin(2C)^ 

- (2(3 - 2^1 + 2C^) cos(2eF) + 8 Cf sin(2^F)) jS f^vfR^- 


(50) 


EFFECT OF ANISOTROPY 


Let us turn to consider the RKKY interaction of the anisotropic Dirac SMs. The effective Hamiltonian for anisotropic 
Dirac SMs in the vicinity of the Weyl node of chirality x takes the form 


^o\ = X (vxkxfTx + Vykyay + vjk^a^ 


(51) 


where Va is the Fermi velocity in the a-direction with a = x,y, z, which generally satisfies the relation It;^; | ^ |uj, | ^ It;^ |. 
The superscript A stands for the anisotropy of Weyl node. To simplify following discussion, we set Va > 0. The 
wave vector ka denotes the derivation from anisotropic Weyl node x- The energy dispersion of Eq. (51) is given 

as = ±\/v'^k'^ + Vyk‘^ + vlk"^. The Green function in momentum space is G“^(k;e) = (e -I- iri)ao — It is 


'fc V ' “v'^v 
instructive to introduce the following associated transformation as 


'^aka ^ VFRa^ — FQa^Fka ^ VaRa^ t 


(52) 
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which leads to two consequences. First, it relates the Green function of an anisotropic Weyl node to that of the 
isotropic one as 


G;^(±R; e) 


d^k 

(27r)3 


Gx(k; e) exp ^ztik • = Aa 


d^k 

(27r)3 


G^(k; e) exp (±ik ■ R) 


AaG(±R; e), 


and 


g (±R,e) 


G+(±R,e) 0 

0 G-(±R,e) 


= Aa 


G+(±R,e) 0 

0 G_(±R,e) 


= XaG (±R,e) 


(53) 


where Aa = Vp/(vxVyVz) measures the anisotropy. Second, it converts a Fermi elliptic sphere to a Fermi sphere but 
preserves its volume, namely, 






We choose R to be aligned to the j-direction {ep = ej) and have 


R = = 


-Rej = 


-R. 


Vp vf 

Substituting this result into the original expression of the RKKY interaction leads to 


HnKKviR) = -tim [ deTr 


{Jtq + Xtx) Si ■ aQ (^R, {Jtq + Xt^) 82 - (tG ^-R,e^ 
A A H'RKKYiR) = Aa Hrkky( — R) , 

\ Vi / 


(54) 


which converts the expression of RKKY interaction of anisotropic Dirac SMs to the counterpart for isotropic Dirac 
SMs. Compared with the isotropic case, the anisotropy not only rescales the amplitude of RKKY interaction by A^, 
but also the effective spatial separation between two impurities by vplvj. Note that this procedure can also be applied 
to the anisotropic Weyl SMs. 


BASIC INTEGRALS 


In this section, we evaluate the following basic integrals 

rip 


In{iF,a)=j' C exjp{ia^)d^, 
J —00 


n = 0 , 1 , 2 . 


We at first consider the integral 


exp(ja^)d^ 

-00 


(55) 


(56) 


which does not converge if a is real. To carry out this integral, we make use of the trick in Ref. [44] by adding a 
negative infinitesimal imaginary part to a, i.e., a ^ a — irj with 77 > 0 , 


T fc rv exp(iaO exp(r 7 ^) 

Ra[^F,a) / exp[z(a - ir])^\d^ = 


i{a — irf) 


-iexp(i<;^)exp(?7^i;-) 


a — IT] 


After taking the limit 77 —>■ O’*", we have 

, , — i exp(m^F) exp( 77 ^F) —iexpUa^p) sin(a^F) cosfa^p) 

Roi^F,a)= lim -^-=-=- i -. 

t;->o+ a — IT] a a a 

It is straightforward to obtain the other basic integrals 

.d cos^a^p) + a^F sm{a^p) .sinia^p) - a^p cosia^p) 

M[^F,a) = -i—Io{^F,a) = - 5 - \-i - r-, 

oa 


(57) 


(58) 


(59) 




d 


R2i^F,a) = (-*) ^ 


da 


{a'^^F — 2) sin(a^F) + 2a^F cos(a^F) . (2 — a'^^p) cos{a^F) + 2a^F sin(a^F) 


(60) 



















